Current-driven switching of magnetisation- theory and experiment 

D. M. Edwards 

Department of Mathematics, Imperial College, London SW7 2BZ, U.K. 

F. Federici 

NEST-INFM and Classe di Scienze, Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy 



I. INTRODUCTION 

Recently there has been a lot of interest in magnetic nanopillars of 10-100 nm in diameter. The pillar is a 
metallic layered structure with two ferromagnetic layers, usually of cobalt, separated by a non-magnetic 
spacer layer, normally of copper. Non-magnetic leads are attached to the magnetic layers so that an 
electric current may be passed through the structure. In the simplest case the pillar may exist in two 
states, with the magnetisation of the two magnetic layers parallel or anti-parallel. The state of a pillar 
can be read by measuring its resistance, this being smaller in the parallel state than in the anti-parallel 
one. This dependence of the resistance on magnetic configuration is the giant magnetoresistance (GMR) 
effect 0. A dense array of these nanopillars could form a magnetic memory for a computer. Normally 
one of the magnetic layers in a pillar is relatively thick and its magnetisation direction is fixed. In order 
to write into the memory the magnetisation direction of the second thinner layer must be switched. This 
might be achieved by a local magnetic field of suitable strength and methods have been proposed |2( for 
providing such a local field by currents in a criss-cross array of conducting strips. However an alternative, 
and potentially more efficient method, proposed by Slonczewski Q makes use of a current passing up 
the pillar itself. Slonczewski's effect relies on "spin transfer" and not on the magnetic field produced by 
the current which in the nanopillar geometry is ineffective. The idea of spin-transfer is as follows. In a 
ferromagnet there are more electrons of one spin orientation than of the other so that current passing 
through the thick magnetic layer (the polarising magnet) becomes spin-polarised. In general its state 
of spin-polarisation changes as it passes through the second (switching) magnet so that spin angular 
momentum is transferred to the switching magnet. This transfer of spin angular momentum is called 
spin-transfer torque and, if the current exceeds a critical value, it may be sufficient to switch the direction 
of magnetisation of the switching magnet. This is called current-induced switching. 

In the next section we show how to calculate the spin-transfer torque for a simple model. 

II. SPIN-TRANSFER TORQUE IN A SIMPLE MODEL 

For simplicity we consider a structure of the type shown in Fig. ^ where p and m are unit vectors in 
the direction of the magnetisations. This models the layered structure of the pillars used in experiments 
but the atomic planes shown are considered to be unbounded instead of having the finite cross-section 
of the pillar. This means that there is translational symmetry in the x and z directions. The structure 
consists of a thick (semi-infinite) left magnetic layer (polarising magnet) , a non-magnetic metallic spacer 
layer, a thin second magnet (switching magnet) and a semi-infinite non-magnetic lead. In the simplest 
model we assume the atoms form a simple cubic lattice, with lattice constant a, and we adopt a one-band 
tight-binding model with hopping Hamiltonian 

H ° = t J2J2 C l ll na C k ll n-la + h.C. (1) 
k|| a n 

Here c£ na creates an electron on plane n with two-dimensional wave- vector k|| and spin er, and t is 
the nearest-neighbour hopping integral. 

In the tight-binding description the operator for spin angular momentum current between planes n—1 
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and n, which we require to calculate spin-transfer torque, is given by 

= - f S (4 u ,Ti,«rT' C * U) n,J " (cfcB.n-LT^Aiy.n-u) 1 + h " C " ( 2 ) 
k ll 

Here er = (cr x , tij,, <t z ) where the components are Pauli matrices. Eq. @ yields the charge current 
jn-i if \ <J is replaced by a unit matrix multiplied by the number e/H, where e is the electronic charge 
(negative). All currents flow in the y direction, perpendicular to the layers, and the components of the 
vector j correspond to transport of x, y and z components of spin. The justification of Eq. J5J) for 
relies on an equation of continuity, as pointed out in Section llVl 

To define the present model completely we must supplement the hopping Hamiltonian Hq by specifying 
the on-site potentials in the various layers. For simplicity we assume the on-site potential for both spins 
in non-magnetic layers, and for majority spin in ferromagnetic layers, is zero. We assume an infinite 
exchange splitting in the ferromagnets so that the minority spin potential in these layers is infinite. Thus 
minority spin electrons are completely excluded from the ferromagnets. Clearly the definition of majority 
and minority spin relate to spin quantisation in the direction of the local magnetisation. We take a = 0, 
so that the magnetisation of the switching magnet is in the z direction and take 6 = ip, where ip is the 
angle between the magnetisations. 

To describe spin transport in the structure we adopt the generalised Landauer approach of Waintal et. 
al. Q. Thus the structure is placed between two reservoirs, one on the left and one on the right, with 
electron distributions characterised by Fermi functions f(u — /il), f(ui — /Ur) respectively. The system is 
then subject to a bias voltage Vy, given by eVt, = /iL — Mr, the difference between the chemical potentials. 
We discuss the ballistic limit where scattering occurs only at interfaces, the effect of impurities being 
negligible. We label the atomic planes so that n = corresponds to the last atomic plane of the polarising 
magnet, the planes of the spacer layer correspond to n = 1,2 • • N and n = N + 1 is the first plane of the 
switching magnet. 

Consider first an electron incident from the left with wave- function \k, maj), where k > 0, which 
corresponds to a Bloch wave |fc) = e lfcT1Q |k|| tt.) with majority spin in the polarising magnet. In this 
notation the label k|| is suppressed. The particle is partially reflected by the structure and finally emerges 
as a partially transmitted wave in the lead, with spin ] corresponding to majority spin in the switching 
magnet. Thus the wave-function is of the form 

\P k ) = \k, maj) + B\ - k, maj) (3) 

in the polariser and 

\L k )=F\k,l) (4) 

in the lead. A majority spin in cither ferromagnet enters or leaves the spacer without scattering, since in 
our simple model there is no potential step. Also the minority spin wave-function entering a ferromagnet 
is zero. The spacer wave-function may therefore be written in two ways: 

\S k ) = F\k, T) + E {c-^ N+1 ^\k, I) - e^+^l - k, I)) (5) 

or 

\Sk) = \k,maj) + B\ — k,maj) + D (\k,min) — | — k,min}) 

= cos 0/2) |fc, t) + sin (V/2) \k, i)+B [cos (V>/2) | - k, f) + sin ty/2) | - k, j>] (6) 
+D [- sin (V-/2) \k, T) + cos 0/2) \k, j) + sin 0/2) | - k, T) - cos 0/2) | - k, j)] . 

On equating coefficients of \k, |), \k, J,), | — k, |), | — k, J.) in expressions © and ® we have four equations 
which may be solved for B, D, E, F. In particular the transmission coefficient T is given by 

T=|F|2= 4cos 2 Q/2)sin 2 fc(iV + l) a 

sin 4 0/2) +4cos 2 0/2)sin 2 fc(A^+ l)a' (> 

Similarly an electron incident from the right with wave- function | — k T) in the lead is partially reflected 
and finally emerges as a partially transmitted wave F'\ — k,maj) in the polarising magnet. It is found 
that F' = F so that the transmission coefficient is the same for particles from left or right. 
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The spin angular momentum current in a particular layer, which we shall denote by S although it need 
not be the spacer layer, is the sum of currents carried by left and right moving electrons. Thus we have 
a Landauer-type formula [f| 

J S = ^E{/ te[(Sk\jn-i\Sk)f(u>-» L ) + (S-k\h-i\S-k)f(u>-»R)}\ (8) 
K ll 

where \Sk), \S-k) are wave-functions in the layer considered corresponding to electrons incident from left 
and right, respectively. Here u>, the energy of the Bloch wave k, is given by the tight-binding formula 

w = Ukii + 2icosA:a (9) 

where Uk» = 2t(cos k x a + cos k z a). We take t < so that positive k corresponds to positive velocity 
h dcj/dk as we have assumed. The current j s in layer S calculated by Eq. JHJ does not depend on the 
particular planes n—l,n between which it is calculated. On changing the integration variable in Eq. (JHJ 
we find 

j- = T E / duJ t J +/( w - Ml) + J-/(w - mr)] (10) 

where 

J± _ (S±k\jn~l\S± k ) 

— 2t sin ka 

Here k = k(ui, k||) is the positive root of Eq. Eq. (f 1 Cjf> may be written as 

j s = f £ / du; {(J+ + J_) [f(cj - ml) + /(w - mr)] + (J+ - J_) [/(w - ml) - /(« - Mr)]} ■ (12) 
* k ii 

Before discussing this spin current we briefly consider the charge current j c , and we denote the analogues 
of J± by JJ. Since the charge current is conserved throughout the structure J+ and JZ can be calculated 
in different ways, e.g. in the lead for Ji and in the polariser for Jf.. Since T = |i^| 2 = |F'| 2 we find 
+ = and for small bias eVb = ml — MR t ne charge current is given by 

f-^-£T (13) 

k ll 

where the transmission coefficient T is given by Eq. (JJJ with fc = k(fi, kii), m being the common chemical 
potential as Vb — + 0. This is the well-known Landauer formula [5j. 
The spin transfer torque on the switching magnet is given by 

T S ' = (jspaccr) — (jlcad), (14) 

where (j spa cer) and (ji ca d) are spin currents in the spacer and lead respectively. For zero bias (ml = 
Mr) there is clearly no charge current in the structure and straight-forward calculation shows that all 
components of spin current in the spacer and the lead vanish, except for a non-zero y-spin current in 
the spacer. There is therefore a non-zero y component of spin-transfer torque acting on the switching 
magnet for zero bias, and its dependence on the angle ip between the magnetisations is found to be 
approximately sin-0. This torque is due to exchange coupling, analogous to an RKKY coupling, between 
the two magnetic layers. This coupling oscillates as a function of spacer thickness and tends to zero as 
the thickness tends to infinity. For finite bias Vb the second term in the integrand of Eq. (|12|l comes 
into play. In general this leads to finite x and y components of T s_t proportional to Vb (for small 14) 
whereas T| _t = 0. However for the special model considered here with infinite exchange splitting in 
both ferromagnets it turns out that TJp* = 0. For this model the only non-zero component of T s_t 
proportional to Vb is found to be 

T B-t = ^f tan ^ (15) 

2e 2 v ; 
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FIG. 1: Schematic picture of a magnetic layer structure for current-induced switching (magnetic layers are darker, 
non-magnetic layers lighter). 

where j c is the charge current given by Eq. I|13|) . 

Slonczewski 3] originally obtained this result for the analogous parabolic band model. From Eqs. 
(|15J) . H13J) and (0) it follows that T| _t contains an important factor sinijj although this does not represent 
the whole angle dependence. Clearly, from Eq. I|15|l . the torque proportional to bias remains finite for 
arbitrarily large spacer thickness, in the ballistic limit. For this model, with infinite exchange splitting, 
the torque is independent of the thickness of the switching magnet. 

From the results of this simple model we can infer a general form of the spin-transfer torque T s_t 
which is independent of the choice of coordinate axes. Thus we write 

T s -*=Tx+T|| (16) 

where 

Tx = (f + Sie7 b )(raxp) 

T|| = g\\eV h m x (p x m). (17) 

With the choice of axes in Fig. fflT|| corresponds to the x component of torque, that is the component 
parallel to the plane containing the magnetisation directions m and p. Similarly Tj_ corresponds to the y 
component of torque, this being perpendicular to the plane of m and p. The modulus of both the vectors 
m x p and m x (p x m) is sin?/;, so that the factors g cx , g± and g\\ are functions of ip which contain 
deviations from the simple sin?/; behaviour. The bias-independent term g ex corresponds to the interlayer 
exchange coupling, as discussed above, and henceforth we assume that the spacer is thick enough for this 
term to be negligible. Sometimes the sinip factor accounts for most of the angular dependence of T± 
and T|| so that g± and g\\ may be regarded as constant parameters for the given structure. In the next 
section we use Eqs. Ijl7|l for the spin-transfer torque in a phenomenological theory of current-induced 
switching of magnetisation. This phenomenological treatment enables us to understand most of the 
available experimental data. It is more usual in experimental works to relate spin-transfer torque to 
current rather than bias. However in theoretical work, based on the Landauer or Keldysh approach, bias 
is more natural. In practice the resistance of the system considered is rather constant (the GMR ratio is 
only a few percent) so that bias and current are in a constant ratio. 

III. PHENOMENOLOGICAL TREATMENT OF CURRENT-INDUCED SWITCHING OF 

MAGNETISATION 

In this section we explore the consequences of the spin-transfer torque acting on a switching magnet 
using a phenomenological Landau Lifshitz equation with Gilbert damping (LLG equation). This is 
essentially a generalisation of the approach used originally by Slonczewski j^] and Sun |(J . We assume 
that there is a polarising magnet whose magnetisation is pinned in the rrz-plane in the direction of a 
unit vector p, which is at general fixed angle 9 to the z-axis as shown in Fig. ^ The pinning of the 
magnetisation of the polarising magnet can be due to its large coercivity (thick magnet) or a strong 
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FIG. 2: Schematic picture of a junction in which current-induced switching is studied experimentally. 

uniaxial anisotropy. The role of the polarising magnet is to produce a stream of spin-polarised electrons, 
i.e. spin current, that is going to exert a torque on the magnetisation of the switching magnet whose 
magnetisation lies in the general direction of a unit vector m. The orientation of the vector m is defined 
by the polar angles a, <f> shown in Fig. ^ There is a non-magnetic metallic layer inserted between the 
two magnets whose role is merely to separate magnetically the two magnetic layers and allow a strong 
charge current to pass. The total thickness of the whole trilayer sandwiched between two non-magnetic 
leads must be smaller than the spin diffusion length l s f so that there are no spin flips due to impurities or 
spin-orbit coupling. A typical junction in which current- induced switching is studied experimentally Q 
is shown schematically in Fig. [21 The thickness of the polarising magnet is 40nm, that of the switching 
magnet 2.5nm and the non-magnetic spacer is 6nm thick. The materials for the two magnets and the 
spacer are cobalt and copper, respectively, which are those most commonly used. The junction cross 
section is oval-shaped with dimensions 60nmxl30nm. A small diameter is necessary so that the torque 
due to the Oersted field generated by a charge current of 10 7 -10 8 A/cm 2 , required for current-induced 
switching, is much smaller than the spin-transfer torque we are interested in. 

The aim of most experiments is to determine the orientation of the switching magnet moment as a 
function of the current (applied bias) in the junction. Sudden jumps of the magnetisation direction, i.e. 
current-induced switching, are of particular interest. The orientation of the switching magnet moment m 
relative to that of the polarising magnet p, which is fixed, is determined by measuring the resistance of 
the junction. Because of the GMR effect, the resistance of the junction is higher when the magnetisations 
of the two magnets are anti-parallel than when they are parallel, In other words, what is observed are 
hysteresis loops of resistance versus current. A typical experimental hysteresis loop of this type [8| is 
reproduced in Fig. [3J It can be seen from Fig. [3]that, for any given current, the switching magnet moment 
is stationary (the junction resistance has a well defined value), i.e. the system is in a steady state. This 
holds everywhere on the hysteresis loop except for the two discontinuities where current-induced switching 
occurs. As indicated by the arrows jumps from the parallel (P) to anti-parallel (AP) configurations of the 
magnetisation, and from AP to P configurations, occur at different currents. It follows that in order to 
interpret experiments which exhibit such hysteresis behaviour, the first task of the theory is to determine 
from the LLG equation all the possible states and then investigate their dynamical stability. At the point 
of instability the system seeks out a new steady state, i.e. a discontinuous transition to a new steady 
state with the switched magnetisation occurs. We have tacitly assumed that there is always a steady 
state available for the system to jump to. There is now experimental evidence that this is not always the 
case. In the absence of any stable steady state the switching magnet moment remains permanently in 
the time-dependent state. This interesting case is implicit in the phenomenological LLG treatment and 
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FIG. 3: Resistance vs current hysteresis loop (after Grollier et al. 



we shall discuss it in detail later. 

In describing the switching magnet by a unique unit vector m, we assume that it remains uniformly 
magnetised during the switching process. This is only strictly true when the exchange stiffness of the 
switching magnet is infinitely large. It is generally a good approximation as long as the switching 
magnet is small enough to remain single domain^so that the switching occurs purely by rotation of the 
magnetisation as in the Stoner-Wohlfarth theory [9j of field switching. This seems to be the case in many 
experiments HHIli [13 

Before we can apply the LLG equation to study the time evolution of the unit vector m in the direction 
of the magnetisation of the switching magnet, we need to determine all the contributions to the torque 
acting on the switching magnet. Firstly, there is the spin-transfer torque T s_t which we discussed in 
Section^ Secondly, there is a torque due to the uniaxial in-plane and easy plane (shape) anisotropies. 
The easy-plane shape anisotropy torque arises because the switching magnet is a thin layer typically 
only a few nanometers thick. The in-plane uniaxial anisotropy is usually also a shape anisotropy arising 
from an elongated cross section of the switching magnet 0. We take the uniaxial anisotropy axis of 
the switching magnet to be parallel to the z-axis of the coordinate system shown in Fig. ^ Since the 
switching magnet lies in the xz-plane, we can write the total anisotropy field as 

H A = H u + H p (18) 

where H u and H p are given by 

B. u = H u0 {m-e z )e z , (19) 



Hp = -if p0 (m • e y )e y . (20) 

Here e x , e y and e z are unit vectors in the directions of the axes shown in Fig. 1. If we write the energy 
of the switching magnet in the anisotropy field as —Ha • (Stot), where (S to t) is the total spin angular 
momentum of the switching magnet, them H u q,H p q which measure the strengths of the uniaxial and 
easy-plane anisotropies have dimensions of frequency. These quantities may be converted to a field in 
tesla by multiplying by h/2fi B — 5.69 x 10~ 12 . 

We are now ready to study the time evolution of the unit vector m in the direction of the switching 
magnet moment. The LLG equation takes the usual form 

dm dm „ . „, 

- +7 mx-.r (21) 

where the reduced total torque T acting on the switching magnet is given by 

r = [- (H A + H ext ) x (S tot ) + Tl + T,|] / |(S tot )| . (22) 



7 



Here H cxt is an external field, in the same frequency units as Ha, and 7 is the Gilbert damping parameter. 
Following Sun [||, Eq. I|21[l may be written more conveniently as 

(l + 7 2 )~r- 7 mxl\ (23) 

It is also useful to measure the strengths of all the torques in units of the strength of the uniaxial 
anisotropy |6J. We shall, therefore, write the total reduced torque T in the form 

T = H uo {(m • e z )m x e z - h p (m ■ e y )m x e y + v^(tp)m x (p x m) + [v±(ip) + h cxt ] m x p} (24) 

where the relative strength of the easy plane anisotropy h p — H p q/H u q and vi\(if)) — vgu^ip), 
v±(ip) = vg±(ip) measure the strengths of the torques T|| and Tj_. The reduced bias is defined by 
v — eVb/(|(S to t)|-ffuo) and has the opposite sign from the bias voltage since e is negative. Thus positive 
v implies a flow of electrons from the polarising to the switching magnet. The last contribution to the 
torque in Eq. (|24|l is due to the external field H cxt with /i cxt = H cxt / H u q. The external field is taken in 
the direction of the magnetisation of the polarising magnet, as is the case in most experimental situations. 

It follows from Eq. (|21(l that in a steady state T = 0. We shall first consider some cases of experimental 
importance where the steady state solutions are trivial and the important physics is concerned entirely 
with their stability. To discuss stability, we linearise Eq. i|23|) . using Eq. Q24|). about a steady state 
solution m = m . Thus 

m = m + £e« + 7^, (25) 

where e a , are unit vectors in the direction m moves when a and 4> are increased independently. The 
linearised equation may be written in the form 

^L = A£ + B m p-=C£, + D v . (26) 
dr dr 

Following Sun we have introduced the natural dimensionless time variable r = i7J u o/(l + 7 2 )- The 
conditions for the steady state to be stable are 

F = A + D sc 0; G = AD - BC ^ (27) 

excluding F = G = 12]. For simplicity we give these conditions explicitly only for the case where 
either Wj^o) = v' ± (ipo) = 0, with -0o = cos _1 (p • mo), or mo = ±p. The case mo = ±p is very common 
experimentally as is discussed below. The stability condition G ^ may be written 

Q 2 v\ + (Qh + cos2a )(Qh + cos 2 a ) + h p {Qh(l - 3 sin 2 <j) sin 2 a ) + cos2a (l - 2 sin 2 a sin 2 </> )} - 
/i 2 sin 2 a sin 2 </> (l - 2 sin 2 O sin 2 a ) ^ 0, (28) 
where vu ~ U||( - 0o)j h = I'xCV'o) + ^oxt and Q = cos'i/'o- The condition F ^ takes the form 

-2(u|| 4- 7/i)Q - 7(cos 2a + cos 2 a ) - "fh p (l - 3 sin 2 4> sin 2 a ) < 0. (29) 

We now discuss several interesting examples, the first of these relating to experiments of Grollier et al. 
|lf)j ] and others. In these experiments the magnetisation of the polarising magnet, the uniaxial anisotropy 
axis and the external field are all collinear (along the in-plane z-axis in our convention) . In this case the 
equation T = 0, with T given by Eq. H24fl . shows immediately that possible steady states are given by 
mo = ±p(«o — 0, 7r), corresponding to the switching magnet moment along the z-axis. These are the 
only solutions when h p = 0. For h p =/= other steady-state solutions may exist but in the parameter 
regime which has been investigated they are always unstable 0|- We shall assume this is always the 
case and concentrate on the solutions mo = ±p. In the state of parallel magnetisation (P) mo = p we 
have Vu = u(7||(0), h = vg±(0) + h cxt , ao = and Q — 1. The stability conditions l(2"5|) and (|2~9"|> become 

[.9H(0)] V + (vg x (0) + h cxt + l) 2 + h p [vgj_(p) + h cxt + 1}>0 (30) 



S-||(0)u + 7 



vg ± (0) + h cxt + 1 + i/ip 



^0. (31) 
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In the state of anti-parallel magnetisation (AP) mo = — p we have = vg^fjr), h — vg±(ir) + h cxt , 
= ir and Q = — 1. The stability conditions for the AP state are thus 

[ 5 ||(7r)] V + (-vg ± (%) - h cxt + l) 2 + h p [-vg±{ir) - h cxt + 1] > (32) 



3||(7I> + 7 



vg±(n) + h cxt -1- -hp 



^ 0. 



(33) 



In the regime of low external field (h cxt ~ 1, i.e. H cxt ~ H u o) we have H p >> H ext (h p fss 100). Eqs. 
301) and lIML'l) may be then approximated by 



vg ± (0) + h cxt + 1 > 



(34) 



ugx(7r) + /i ext -l<0. (35) 

Equation l)34|) corresponds to P stability and l|35ll to AP stability. It is convenient to define scalar 
quantities T±, Th by T± = g\_{ip) sin ip, Tji — g\\(ip) sintp, these being scalar components of spin-transfer 
torque in units of eVu fcf. Eq. JEJl). Then gi(0) = [dT i /dV']^=o and gi(n) = — [dTj/d^]^, =7r with i =_L, ||. 
Model calculations 13] show that both g± and g\\ can be of either sign, although positive values are more 
common. Also there is no general rule about the relative magnitude of gi(0) and gi(ir). 

We now illustrate the consequences of the above stability conditions by considering two limiting cases. 
We first consider the case g±(ip) = 0, gu > 0, as assumed by Grollier et. al 10] in the analysis of their data. 
In Fig. 4 we plot the regions of P and AP stability deduced from Eqs. l|^T|l . (|3*3^l - ll^|) . in the (u,ft, ext )-plane. 
Grollier et al. plot current instead of bias but this should not change the form of the figure. Theirs is 
rather more complicated, owing to a less transparent stability analysis with unnecessary approximation. 
The only approximations made above, to obtain Eqs. H34JI and lj35(l . can easily be removed, which results 
in the critical field lines h cxt = ±1 acquiring a very slight curvature given by h ex t «1+ [v g» (ir)] 2 /h p and 
h^t ~ — 1 — [vgii (0)] 2 /h p . The critical biases in the figure are give by 



^AP^P 



1 + ~ ^cxt 



Vp^AP = -7 



1 

2 >Lp 



1 + -h n + h ext 



/siiW 

/.9||(0). (36) 



A downward slope from left to right of the corresponding lines in Fig. 0|is not shown there. Since the 
damping parameter 7 is small (7 ~ 0.01) this downward slope of the critical bias lines is also small. From 
Fig. 01 we can deduce the behaviour of resistance versus bias in the external field regimes \h cxt \ < 1 and 
|/lext| > 1. 

Consider first the case |/i e xt| < 1- Suppose we start in the AP state with a bias v = which is gradually 
increased to vap^p- At this point the AP state becomes unstable and the system switches to the P state 
as v increases further. On reducing v the hysteresis loop is completed via a switch back to the AP state 
at the negative bias vp^ap- The hysteresis loop is shown in Fig. Eta)- The increase in resistance R 
between the P and AP states is the same as would be produced by varying the applied field in a GMR 
experiment. Now consider the case h cx t < — 1- Starting again in the AP state at v — we see from Fig. 0] 
that, on increasing v to wap^p, the AP state becomes unstable but there is no stable P state to switch to. 
This point is marked by an asterisk in Fig. Etb). For v > vap^p, the moment of the switching magnet 
is in a persistently time-dependent state. However, if v is now decreased below up^ap the system homes 
in on the stable AP state and the overall behaviour is reversible, i.e. no switching and no hysteresis 
occur. When h cxt > 1 similar behaviour, now involving the P state, occurs at negative bias, as shown in 
Fig. Efb). The dashed curves in Fig. E^b) show a hypothetical time-averaged resistance in the regions of 
time-dependent magnetisation. As discussed later time-resolved measurements of resistance suggest that 
several different types of dynamics can occur in these regions. 

It is clear from Fig. EJa) that the jump AP— *P always occurs for positive bias v, which corresponds to 
flow of electrons from the polarising to the switching magnet. This result depends on the assumption that 
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FIG. 5: (a) Hysteresis loop of resistance vs bias for |/i cx t| < 1; (b) Reversible behaviour (no hysteresis) for 
|/iext| < — 1 (upper curve) and /i ex t > 1 (lower curve). The dashed lines represent hypothetical behaviour of 
average resistance in regions of Fig. [I] marked "both unstable" where no steady states exist. 



9|| > 0; if gn < it is easy to see that the sense of the hysteresis loop is reversed and the jump P-^AP 
occurs for positive v. To our knowledge this reverse jump has never been observed, although < 
can occur in principle and is predicted theoretically for the Co/Cu/Co(lll) system with a switching 
magnet consisting of a single atomic plane of Co. It follows from Eq. (|36|) that |vp->Ap/f ap-^p| — 
l5||( 7r )/ff||(0)l m zero external field. Experimentally this ratio, essentially the same as the ratio of critical 
currents, may be considerably less than 1 (e.g. < 0.5 0), greater than 1 (e.g. « 2 or close to 1 
Usually the field dependence of the critical current is found to be stronger than that predicted by Eq. 

We now discuss the reversible behaviour shown in Fig. GJb) which occurs for \h cxt \ > 1. The transition 
from hysteretic to reversible behaviour at a critical external field seems to have been first seen in pillar 
structures by Katine et al. |l4j . Curves similar to the lower one in Fig. EJb) are reported with |up_>ap| 
increasing with increasing /i cx t, as expected from Eq. (|36|l . Plots of the differential resistance dV/dl show 
a peak near the point of maximum gradient of the dashed curve. Similar behaviour has been reported 
by several groups [l5l IT^ . Il7| . It is particularly clear in the work of Kiselev at al. 0] that the transition 
from hysteretic behaviour (as in Fig. Efa)) to reversible behaviour with peaks in dV/dl occurs at the 
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coercive field 600 Oe of the switching layer (h cxt = 1)- The important point about the peaks in dV/dl 
is that for a given sign of /i cxt they only occur for one sign of the bias. This clearly shows that this 
effect is due to spin-transfer and not to Oersted fields. Myers et al. show a current-field stability 
diagram similar to the bias-field one of Fig. 0]with a critical field of 1500 Oe. They examine the time 
dependence of the resistance at room temperature with the field and current adjusted so that the system 
is in the "both unstable" region in the fourth quadrant of Fig. 01 but very close to its top left-hand 
corner. They observe telegraph- noise- type switching between approximately P and AP states with slow 
switching times in the range 0.1-10 s. Similar telegraph noise with faster switching times was observed by 
Urazhdin et a l. [lq at current and field close to a peak in dV/dl. In the region of P and AP instability 
Kiselev et al. |l5l| and Pufall et al. report various types of dynamics of p recessional type and random 
telegraph switching type in the microwave Ghz regime. Kiselev et al. 15] propose that systems of the 
sort considered here might serve as nanoscale microwave sources or oscillators, tunable by current and 
field over a wide frequency range. 

We now return to the stability conditions (|31J) . (|33|I - (|35JI and consider the case of g±(i[>) ^ but h CKt = 0. 
These conditions of stability of the P state may be written approximately, remembering that 7 << 1, 
h p >> 1, as 

vg±(0)>-l, W3||(0)>-^7V (37) 
The conditions for stability of the AP state are 

vgx(ir)<l, vg { \ (tt) < ~lK- ( 38 ) 

In Fig. [S] we plot the regions of P and AP stability, assuming <7_l(0) = g±(Tf) — g± and g» (0) = gn (tt) = gu 
for simplicity. We also put r = g±/g\\. For r > we find the normal hysteresis loop as in Fig. Efa) 
if we plot R against vg\\ (valid for either sign of <7||). In Fig. we plot the hysteresis loops for the 
cases r c < r < and r < r c , where r c = —2/('yhp) is the value of r at the point X in Fig. The 
points labelled by asterisks have the same significance as in Fig. OJb). If in Fig. |7{a) we increase vg\\ 
beyond its value indicated by the right-hand asterisk we move into the "both-unstable" region where the 
magnetisation direction of the switching magnet is perpetually in a time-dependent state. Thus negative 
r introduces behaviour in zero applied field which is similar to that found when the applied field exceeds 
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FIG. 7: Hysteresis loop for (a) r c < r < 0; (b) r < r c 



the coercive field of the switching magnet for r — 0. This behaviour was predicted by Edwards et al. 
|l3j. in particular for a Co/Cu/Co(lll) system with the switching magnet consisting of a Co monolayer. 
Zimmler et al. use methods similar to the ones described here to analyse their data on a Co/Cu/Co 
nanopillar and deduce that g\\ > 0, r — g±/g\\ ~ —0.2. It would be interesting to carry out time-resolved 
resistance measurements on this system at large current density (corresponding to vg± < — 1) and zero 
external field. 

So far we have considered the low-field regime (-ff cx t ~ coercive field of switching magnet) with both 
magnetisations and the external field in-plane. There is another class of experiments in which a high 
field, greater than the demagnetising field (> 2T), is applied perpendicular to the plane of the layers. The 
magnetisation of the polarising magnet is then also perpendicular to the plane. This is the situation in 
the early experiments where a point contact was employed to inject high current densities into magnetic 
multilayers [201 l2ll l22| . In this high- field regime a peak in the differential resistance dV/dl at a critical 
current was interpreted as the onset of current-induced excitation of spin waves in which the spin-transfer 
torque leads to uniform precession of the magnetisation [HEH, HI. No hysteretic magnetisation reversal 
was observed and it seemed that the effect of spin-polarised current on the magnetisation is quite different 
in the low- and high-field regimes. Recently, however, Ozyilmaz et al. [24) have studied Co/Cu/Co 
nanopillars (« lOOnm in diameter) at T = 4.2K for large applied fields perpendicular to the layers. They 
observe hysteretic magnetisation reversal and interpret their results using the Landau-Lifshitz equation. 
We now give a similar discussion within the framework of this section. 

Following Ozyilmaz et at, we neglect the uniaxial anisotropy term in Eq. (|24(l for the reduced torque 
r while retaining H u q as a scalar factor. Hence 



T = H u0 {[/i ex t + v±(ij}) - h p cos V>] m x p + v\\(il))m x (p x m)} 



(39) 



where p is the unit vector perpendicular to the plane. When uii (ip) ^ the only possible steady-state 
solutions of r = are mo = ±p. On linearizing Eq. [23 about mo as before we find that the condition 
G ^ is always satisfied. The second stability condition F < becomes 



[w||(V>o) +l(v±(ip ) + h cxt - hp)] cos -00 > 



(40) 



where ipo — cos 1 (mo • p). Applying this to the P state (^o = 0) and the AP state (ipo ~ it) we obtain 
the conditions 



V > "{{hp - /lext)/#(0) 



(41) 



v < -j(h p + h cxt )/g(Tr), 



(42) 



where the first condition applies to the P stability and the second to the AP stability. Here g(ip) — g\\ + 
r yg±(tp). The corresponding stability diagram is shown in Fig. |S| where we have assumed gin) > g(0) > 
for definiteness. The boundary lines cross at h CKt — h Cl where h c = hp{g(ir) + g(0)]/{g(ir) — g(0)]. This 
analysis is only valid for fields larger than the demagnetising field (h cxt > h p ) and we see from the figure 
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FIG. 8: Bias-field stability diagram for large external field (/i ex t > h p ) perpendicular to the layers. 

that for h cyLt > h c hysteretic switching occurs. This takes place for only one sign of the bias (current) 
and the critical biases (currents) increase linearly with h ext as does the width of the hysteresis loop 
l^p— >ap — fAP->p|- This accords with the observations of Ozyilmaz et al. The critical currents are not 
larger than those in the low- field or zero-field regimes (cf. Eqs. (f 4 1|) . 142JI with Eq. (j2EJ)) and yet the 
magnetisation of the switching magnet can be switched against a very large external field. However, in 
this case the AP state is only stabilised by maintaining the current. 

The experiments on spin transfer discussed above have mainly been carried out at constant temperature, 
typically 4.2K or room temperature. The effect on current-driven switching of varying the temperature 
has recently been studied by several groups 0, 0, |25| . The standard Neel-Brown theory of thermal 
switching [26j does not apply because the Slonczewski in-plane torque is not derivable from an energy 
function. Li and Zhang }27j have generalised the standard stochastic Landau-Lifschitz equation, which 
includes white noise in the effective applied field, to include spin transfer torque. In this way they have 
successfully interpreted some of the experimental data. A full discussion of this work is outside the 
scope of the present review. However it should be pointed out that in addition to the classical effect of 
white noise there is an intrinsic temperature dependence of quantum origin. This arises from the Fermi 
distribution functions which appear in expressions for the spin-transfer torque (see Eqs. (|14|l and l|12|l L 

So far we have discussed steady-state solutions of the LLG equation l|23|) . It is important to study the 
magnetisation dynamics of the switching layer in the situation during the jumps AP— >P and P— >AP of the 
hysteresis curve in zero external field, and secondly under conditions where only time-dependent solutions 
are possible, for example in the regions of sufficiently strong current and external field marked "both 
unstable" in Fig. 21 The first situation has been studied by Sun [(J, assuming single-domain behaviour 
of the switching magnet, and by Miltat et al. |28| with more general micromagnetic configurations. Both 
situations have been considered by Li and Zhang 29]. In the second case they find precessional states, 
and the possibility of "telegraph noise" at room temperature, as seen experimentally in Refs. [TH ITi|. 
Switching times (AP^P and P— »AP) are estimated to be of the order Ins. Micromagnetic simulations 
[28l | indicate that the Oersted field cannot be completely ignored for typical pillars with diameter of the 
order of lOOnm. 

Finally, in this section, we briefly discuss some practical considerations which may ultimately decide 
whether current-induced switching is useful in spintronics. Sharp switching, with nearly rectangular 
hysteresis loops, is obviously desirable and this demands single-domain behaviour. In experiments on 
nanopillars of circular cross section [llj multidomain behaviour was observed with the switching tran- 
sition spread over a range of current. Subsequently the same group found sharp switching in pillars 
whose cross-section was an elongated hexagon, which introduces strong uniaxial in-plane shape anisotropy. 
It was known from earlier magnetisation studies of nanomagnet arrays |3Cj that such a shape anisotropy 
can result in single domain behaviour. A complex switching transition need not necessarily indicate mul- 
tidomain behaviour. It could also arise from a marked departure of Tj_(ip) and/or Tji ("0) from sinusoidal 
behaviour, such as occurs near ip — it in calculations for Co/Cu/Co(lll) with two atomic planes of Co 
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in the switching magnet (see Fig. EJb)). In the calculations of the corresponding hysteresis loops (Fig. 
II ip the torques were approximated by sine curves but an accurate treatment would certainly complicate 
the AP^P transition which occurs at negative bias in Fig. Illf b). Studies of this effect are planned. 

The critical current density for switching is clearly an important parameter. From Eq. (|36|l the critical 
reduced bias for the P^AP transition is to a good approximation given by — jh p /[2gn(0)]. Using the 
definitions of reduced quantities given after Eq. (|24(l , we may write the actual critical bias in volts as 

Vp^ap = Af 7 M s iV[2 5 ||(0)|e|], (43) 

where M is the number of atomic planes in the switching magnet, M s is the average moment (J/T) of 
the switching magnet per atomic plane per unit area, and Hd = HH v o/ (2/xb) is the easy-plane anisotropy 
field in tesla. As expressed earlier g\\ (0) = (dTji/dV>)i/>=o where the torque Ty is per unit area in units of 
eVe- (The calculated torques in Figs. I§1 and ITUl of Sec. 0are per surface atom so that if these are used 
to determine g\\(0) in Eq. (|43|l M s must be taken per surface atom.) 

An obvious way to reduce the critical bias, and hence the critical current, is to reduce M, the thickness 
of the switching magnet. Calculations show 0] (see also Fig. I10|) that gn does not decrease with M and 
may, in fact, increase for small values such as M = 2. Careful design of the device might also increase 
<7||(0) beyond the values (< 0.01 per surface atom) which seem to be obtainable in simple trilayers |l3j |. 
Jiang et al. 0, 112, have studied various structures in which the polarising magnet is pinned by an 
adjacent antiferromagnet (exchange biasing) and in which a thin Ru layer is incorporated between the 
switching layer and the lead. Critical current densities of 2 x 10 6 Acm~ 2 have been obtained which are 
substantially lower than those in Co/Cu/Co trilayers. Such structures can quite easily be investigated 
theoretically by the methods of Section Ivl 

Decreasing the magnetisation M s , and hence the demagnetising field (oc i?d), would be favourable 
but g<i then tends to decrease also A possible way of decreasing without decrea sing local 

magnetic moments in the system is to use a synthetic ferrimagnet as the switching magnet 33]. The 
Gilbert damping factor 7 is another crucial parameter but it is uncertain whether this can be decreased 
significantly. However, the work of Capelle and Gyorffy |34| is an interesting theoretical development. The 
search for structures with critical current densities low enough for use in spintronic devices (10 5 Acm~ 2 
perhaps) |35j | is an enterprise where experiment and quantitative calculations |l3j | should complement 
each other fruitfully. 



IV. QUANTITATIVE THEORY OF SPIN-TRANSFER TORQUE 
A. General principles 

To put the phenomenological treatment of Sec. IIIII on a first-principle quantitative basis we must 
calculate the spin-transfer torques (Eqs. (|TT|) in a steady state for real systems. For this purpose it is 
convenient to describe the magnetic and nonmagnetic layers of Fig. fflby tight-binding models, in general 
multiorbital with s, p, and d orbitals whose one-electron parameters are fitted to first-principle bulk band 
structure [3(|. The hamiltonian is therefore of the form 

H = H Q + H int + ff anis (44) 

where the one-electron hopping term Hq is given by 

H 0=Y1 Yl W,rw( k ||) c k || m/ JCT c k||n^' ( 45 ) 

where c\ mfia creates an electron in a Bloch state, with in-plane wave vector k|| and spin cr, formed 

from a given atomic orbital /1 in plane m. Eq. 1451 generalises the single orbital eq. Q. H[ nt is an 
on-site interaction between electrons in d orbitals which leads to an exchange splitting of the bands in 
the ferromagnets and is neglected in the spacer and lead. Finally, i? an is contains anisotropy fields in the 
switching magnet and is given by 



-ffanis — / y S n • Ha, 

71 



(46) 
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where S„ is the operator of the total spin angular momentum of plane n and Ha is given by Eqs. (|18|) - 
(|2~U)l with the unit vector m in the direction of ^ n {S n ), where (S„) is the thermal average of S„. We 
assume here that the anisotropy fields H u q,H p are uniform throughout the switching magnet but we could 
generalise to include, for example, a surface anisotropy. 

In the tight-binding description, the spin angular momentum operator S„ is given by 

and the corresponding operator for the spin angular momentum current between planes n — 1 and n is 

jn-l = -\ H Y1 *( k l|W-lA< ( C l||n<4'4||r»/i) CT n-lu T , C fc|| T + h.C. , (48) 

k||/*n 

which generalises the single orbital expression @. The rate of change of S„ in the switching magnet is 
given by 

iH—^ = [S n ,H ] + [S n ,H aIlis \. (49) 

This results holds since the spin operator commutes with the interaction hamiltonian Hi nt . 
It is straightforward to show that 

[S„,fr ]=i»(j n _i-j„), (50) 

and 



[S n ,H anis ] = -ih(U A xS n ). (51) 
On taking the thermal average, Eq. I|49|l becomes 

<^> = <j«-i> - On) - Ha x (Stat), (52) 

This corresponds to an equation of continuity, stating that the rate of change of spin angular momentum 
on plane n is equal to the difference between the rate of flow of this quantity onto and off the plane, plus 
the rate of change due to precession around the field Ha- When Eq. (|52|l is summed over all planes in 
the switching magnet we have 

^(S tot )=T s - t -H A x (Stat), (53) 

where the total spin-transfer torque T s_t is given by Eq. (|14|) and (Stot) is the total spin angular 
momentum of the switching magnet. Equation (|53|l is equivalent to Eq. (|21|l . for zero external field, in 
the absence of damping. Equation (|14(l shows how T s-t required for the phenomenological treatment of 
Sec. IIIII is to be determined from the calculated spin currents in the spacer and lead. As discussed in 
Sec. IIIII the magnetization of a single-domain sample is essentially uniform and the spin-transfer torque 
T s_t depends on the angle ip between the magnetisations of the polarising and switching magnets. 

To consider time-dependent solutions of Eq. H21(l it is necessary to calculate T s_t for arbitrary angle 
if) and for this purpose Ha can be neglected. To reduce the calculation of the spin-transfer torque to 
effectively a one-electron problem, we replace H lvA by a selfconsistent exchange field term — S„ • A„, 
where the exchange field A n should be determined selfconsistently in the spirit of an unrestricted Hartree- 
Fock (HF) or local spin density (LSD), approximation. The essential selfconsistency condition in any HF 
or LSD calculation is that the local moment (S„) in a steady state is in the same direction as A„. Thus 
we require 

A n x (S„) = (54) 

for each atomic plane of the switching magnet. It is useful to consider first the situation when there is 
no applied bias and the polarising and switching magnets are separated by a spacer which is so thick 
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that the zero-bias oscillatory exchange coupling [37j is negligible. In that case we have two independent 
magnets and the selfconsistent exchange field in every atomic plane of the switching magnet is parallel 
to its total magnetisation which is uniform and assumed to be along the z-axis. Referring to Fig. ^thc 
selfconsistent solution therefore corresponds to uniform exchange fields in the polarising and switching 
magnets which are at an assumed angle ip = 9 with respect to one another. 

When a bias Vt, is applied , with a uniform exchange field A = Ae z in the switching magnet imposed, 
the calculated local moments (S„) will deviate from the z-direction so that the solution is not selfcon- 
sistent. To prepare a selfconsistent state with A and all (S„) = (S) in the z-direction it is necessary to 
apply fictitious constraining fields H„ of magnitude proportional to Vt>. The local field for plane n is thus 
A + H n but to calculate the spin currents in the spacer and lead, and hence T s_t from Eq. I|14|) . the 
fields H„, of the order of Vb, may be neglected compared with A. Although the fictitious constraining 
fields H„ need therefore never be calculated, it is interesting to see that they are in fact related to T s_t . 
For the constrained self-consistent steady state ((S n ) = (S), (S„) = 0) in the presence of the constraining 
fields, with Ha neglected as discussed above, it follows from Eq. (|52() that 

(jn-l) - On) = (A + Hn) X (S) = H„ X (S), (55) 

where the local field A + H„ replaces Ha- On summing over all atomic planes n in the switching magnet 
we have 

T S -* = (j spaccl .) - (head) = H » X < S >- ( 56 ) 

n 

Thus, as expected, in the prepared state with a given angle ip between the magnetisations of the magnetic 
layers the spin-transfer torque is balanced by the total torque due to the constraining fields. 

In the simple model of Section [n] with infinite exchange splitting in the magnets, the local moment is 
constrained to be in the direction of the exchange field so the question of selfconsistency is not raised. 

The main conclusion of this Section is that the spin-transfer torque for a given angle ip between 
magnetisations may be calculated using uniform exchange fields making the same angle with one another. 
Such calculations are described in Sec. ITTlandlVl The use of this spin-transfer torque in the LLG equation 
of Section ITTT1 completes what we shall call the "standard model" (SM). It underlies the original work 
of Slonczewski 3J and most subsequent work. The spin-transfer torque calculated in this way should 
be appropriate even for time-dependent solutions of the LLG equation. This is based on the reasonable 
assumption that the time for the electronic system to attain a " constrained steady state" with given ip is 
short compared with the time-scale («lns) of the macroscopic motion of the switching magnet moment. 

Although the SM is a satisfactory way of calculating the spin-transfer torque its lack of selfconsistency 
leads to some non-physical concepts. The first of these is the "transverse spin accumulation" in the 
switching magnet |3Jj|39|, This refers to the deviations of local moments (S n ) from the direction of the 
exchange field, assumed uniform in the SM. In a self-consistent treatment such deviations do not occur 
because the exchange field is always in the direction of the local moment. A related non-physical concept 
is the " spin decoherence length" over which the spin accumulation is supposed to decay (H, |3^| , More 
detailed critiques of these concepts are given elsewhere [T3L l4fj| . 



B. Keldysh formalism for fully realistic calculations of the spin-transfer torque 

The wave-function approach to spin-transfer torque described in Section ^ is difficult to apply to 
realistic multiorbital systems. For this purpose Green functions are much more convenient and Keldysh 
|4lj | developed a Green function approach to the non-equilibrium problem of electron transport. In this 
section we apply this method to calculate spin currents in a magnetic layer structure, following Edwards 
et al. 13]. 

The structure we consider is shown schematically in Fig. ^ It consists of a thick (semi-infinite) left 
magnetic layer (polarising magnet), a nonmagnetic metallic spacer layer of N atomic planes, a thin 
switching magnet of M atomic planes, and a semi-infinite lead. The broken line between the atomic 
planes n — 1 and n represents a cleavage plane separating the system into two independent parts so that 
charge carriers cannot move between the two surface planes n — 1 and n. It will be seen that our ability 
to cleave the whole system in this way is essential for the implementation of the Keldysh formalism. This 
can be easily done with a tight-binding parametrisation of the band structure by simply switching off the 
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matrix of hopping integrals t^n-i^ between atomic orbitals v, fi localised in planes n — 1 and n. We 
therefore adopt the tight-binding de scrip tion with the Hamiltonian defined by Eqs. (|44I47() . 

To use the Keldysh formalism |4lL l42ll4^ | to calculate the charge or spin currents flowing between the 
planes n — 1 and n, we consider an initial state at time r = — oo in which the hopping integral t nv ^i—Xn 
between planes n — 1 and n is switched off. Then both sides of the system are in equilibrium but with 
different chemical potentials /^l on the left and /zr on the right, where /il — fJ-R = &Vh- The interplane 
hopping is then turned on adiabatically and the system evolves to a steady state. The cleavage plane, 
across which the hopping is initially switched off, may be taken in either the spacer or in one of the 
magnets or in the lead. In principle, the Keldysh method is valid for arbitrary bias Vt, but here we 
restrict ourselves to small bias corresponding to linear response. This is always reasonable for a metallic 
system. For larger bias, which might occur with a semiconductor or insulator as spacer, electrons would 
be injected into the right part of the system far above the Fermi level and many-body processes neglected 
here would be important. Following Keldysh 0,^3 , we define a two-time matrix 

G+ L (r,r')=i(4(r')c R (r)} ) (57) 

where R = (n,v,<j') and L = (n — 1,/J,, cr), and we suppress the fc|| label. The thermal average in Eq. 

(|57|) is calculated for the steady state of the coupled system. The matrix Gjjt has dimensions 2m x 2m 
where m is the number of orbitals on each atomic site, and is written so that the m x m upper diagonal 
block contains matrix elements between f spin orbitals and the m x m lower diagonal block relates to 
I spin. 2m x 2m hopping matrices £lr and £rl are written similarly and in this case only the diagonal 
blocks are nonzero. If we denote £lr by t, then £rl = t We also generalise the definition of cr so that its 
components are now direct products of the 2x2 Pauli matrices a Xl a y , a z , and the m x m unit matrix. 
The thermal average of the spin current operator, given by Eq. (|4*§1> . may now be expressed as 

(j«-i> = \ E Tr { t G RL (r, r) * - G+ R (r, r)tt] cr} . (58) 

k|| 

Introducing the Fourier transform G + (o>) of G + (t, r') , which is a function of r — r', we have 

<]»-!> = / ^lr{[G+,( W )t-G+ a ( W )tV}. (59) 

The charge current is given by Eq. (|59|l with i<r replaced by the unit matrix multiplied by e/H. 
Following Keldysh |4ll |42| we now write 

GabM = ^ab + G ab -G ab ), (60) 
where the suffices A and B are either R or L. _Fab(^) is the Fourier transform of 

F AB (ry) = -i([ CA (r),4(r')]_) (61) 

and G a , G r are the usual advanced and retarded Green functions (SJ. Note that in an d the 
definitions of G a and G r are interchanged and that in the Green function matrix defined by these authors 
G + and G~ should be interchanged. 

Charge and spin current are related by Eqs. l|5§|l and (ffTTTf) to the quantities G a , G r and Fab- The 
latter are calculated for the coupled system by starting with decoupled left and right systems, each in 
equilibrium, and turning on the hopping between planes L and R as a perturbation. Hence, we express 
G a , G r and Fab hi terms of retarded surface Green functions 5l = 9r = 9rr for the decoupled 
equilibrium system. It is then found [Iflj that the spin current between the planes n — 1 and n can be 
written as the sum (j ra _i) = (j n -i)i + (jn— 1)2, where the two contributions to the spin current (j n )i> 
(jn) 2 are given by 

(jn-i>i = 7-E / ~ A)tr]\f(w - /M h ) + f(w - pa)]. (62) 
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(j„-i>2 = ^ E / dcuWTr {[gtfABgltf - AB + ~(A + S)]<r| [f(u> - ,i L ) - /(w - (63) 
k n 

Here, A = [1 — gp^f/Li] 1 , B = [1 — glf^i] -1 , anc [ as m Section |H] /(w — m) i s the Fermi function 
with chemical potential fi and /il — /^r = eVt,. In the linear-response case of small bias which we are 
considering, the Fermi functions in Eq. 1(63|) are expanded to first order in Vb- Hence the energy integral is 
avoided, being equivalent to multiplying the integrand by eVb and evaluating it at the common zero-bias 
chemical potential fj,Q- 

It can be seen that Eqs. ((62(1 and 1(63(1 . which determine the spin and the charge currents, depend on just 
two quantities, i.e. the surface retarded one-electron Green functions for a system cleaved between two 
neighbouring atomic planes. The surface Green functions can be determined without any approximations 
by the standard adlayer method (see e.g. [sHES]) for a fully realistic band structure. 

We first note that there is a close correspondence between Eqs. 1)62(1. 163(1 and the generalised Landauer 
formula 1(12|) . The first term in Eq. ((12|) corresponds to the zero-bias spin current given by Eq. 

((62(1 . When the cleavage plane is taken in the spacer, the spin current { j n — l ) l determines the oscillatory 
exchange coupling between the two magnets and it is easy to verify that the formula for the exchange 
coupling obtained from Eq. 1(62(1 is equivalent to the formula used in previous total energy calculations of 
this effect |36l l37j . The contribution to the transport spin current given by Eq. ((63(1 clearly corresponds 
to the second term in the Landauer formula 1(12(1 which is proportional to the bias in the linear response 
limit. Placing the cleavage plane first between any two neighbouring atomic planes in the spacer and 
then between any two neighbouring planes in the lead, we obtain from Eq. (JH3 the total spin-transfer 
torque T s ~* of Eq. (JT3J in Section HJ 

The equivalence of the Keldysh and Landauer methods has been demonstrated by calculating the 
currents 1(62(1 and ((63(1 analytically for the simple single orbital model of Section [H] The results of that 
section, such as Eq. (|T3|l are reproduced [l3| . 



V. QUANTITATIVE RESULTS FOR Co/Cu/Co(lll) 



We now discuss the application of the Keldysh formalism to a real system. In particular we consider a 
realistic multiorbital model of fee Co/Cu/Co(lll) with tight-binding parameters fitted to the results of 
the first-principles band structure calculations, as described previously [3fj l37| . 

Referring to Fig. ^ the system considered by Edwards et al. [13| consists of a semi-infinite slab of Co 
(polarising magnet), the spacer of 20 atomic planes of Cu, the switching magnet containing M atomic 
planes of Co, and the lead which is semi-infinite Cu. The spacer thickness of 20 atomic planes of Cu 
was chosen so that the contribution of the oscillatory exchange coupling term is so small that it can be 
neglected. The spin currents in the right lead and in the spacer were determined from Eq. IpjS)) . Figure 
E[a),(b) shows the angular dependences of Tu, T± for the cases M = 1 and M — 2. respectively. For 
the monolayer switching magnet, the torques Tji and T± are equal in magnitude and they have opposite 
sign. However, for M = 2, the torques have the same sign and T± is somewhat smaller than Tji. A 
negative sign of the ratio of the two torque components has important and unexpected consequences for 
hysteresis loops as already discussed in Section IIHI It can be seen that the angular dependence of both 
torque components is dominated by a sin^ factor but distortions from this dependence are clearly visible. 
In particular, the slopes at if) = and ip = it are quite different. As pointed out in Section IlIII this is 
important in the discussion of the stability of steady states and leads to quite different magnitudes of the 
critical biases Vp — > Vap and Vap — » Vp. 

In Fig. E| we reproduce the dependence of T± and Tji on the thickness of the Co switching magnet. 
It can be seen that the out-of-plane torque Tj_ becomes smaller than Tji for thicker switching magnets. 
However, T± is by no means negligible (27% of Tii) even for a typical experimental thickness of the 
switching Co layer of ten atomic planes. It is also interesting that beyond the monolayer thickness, the 
ratio of the two torques is positive with the exception of M — 4. 

The microscopically calculated spin-transfer torques for Co/Cu/Co(lll) were used by Edwards et al. 
|13| as an input into the phenomenological LLG equation. For simplicity the torques as functions ofi/i were 
approximated by sine curves but this is not essential. The LLG equation was first solved numerically to 
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FIG. 9: Dependence of the spin-transfer torque Tj| and T± for Co/Cu/Co(lll) on the angle ip. The torques per 
surface atom are in units of eVh- Figure (a) is for M = 1, and (b) for M = 2 monolayers of Co in the switching 
magnet. 
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FIG. 10: Dependence of the spin-transfer torque Tj| and T± for Co/Cu/Co(lll) on the thickness of the switching 
magnet M for ip — ir/3. The torques are in units of eVh- 



determine all the steady states and then the stability discussion outlined in the phenomenological section 
was applied to determine the critical bias for which instabilities occur. Finally, the ballistic resistance of 
the structure was evaluated from the real-space Kubo formula at every point of the steady state path. 
Such a calculation for the realistic Co/Cu system then gives hysteresis loops of the resistance versus bias 
which can be compared with the observed hysteresis loops. The LLG equation was solved including a 
strong easy-plane anisotropy with h p = 100. If we take H u o = 1.86 x 10 9 sec _1 , corresponding to a uniaxial 
anisotropy field of about 0.01T, this value of h p corresponds to the shape anisotropy for a magnetisation 
of 1.6 x 10 6 A/m, similar to that of Co Also a realistic value Q of the Gilbert damping parameter 
7 = 0.01 was used. Finally, referring to the geometry of Fig. ^ two different values of the angle 9 were 
employed in these calculations: 9 = 2rad and 9 = 3rad, the latter value being close to the value of ir 
which is realised in most experiments. 

We first reproduce in Fig. ^2 the hysteresis loops for the case of Co switching magnet consisting of 
two atomic planes. We note that the ratio r = T±/T\\ w 0.65 deduced from Fig. El is positive in this 
case. Fig. Illf a) shows the hysteresis loop for 9 = 2 and Fig. II If that for 9 = 3. The hysteresis loop 
for 9 = 3 shown in Fig. Illf b) is an illustration of the stability scenario in zero applied field with r > 




FIG. 11: Resistance of the Co/Cu/Co(lll) junction as a function of the applied bias, with M = 2 monolayers of 
Co in the switching magnet, (a) is for 6 = 2rad and (b) for 6 = 3rad. 




FIG. 12: Resistance of the Co/Cu/Co(lll) junction as a function of the applied current, with M = 1 monolayer 
of Co in the switching magnet, (a) is for 6 — 2rad and (b) for 9 = 3rad. 



discussed in Section IIIII As pointed out there the hysteresis curve is that of Fig. [Sfa) which agrees with 
Fig. IllT b) when we remember that the reduced bias used in Fig. 0has the opposite sign from the bias in 
volts used in Fig. ^] It is rather interesting that the critical bias for switching is « 0.2mV both for 9 = 2 
and 9 = 3. When this bias is converted to the current density using the calculated ballistic resistance of 
the junction, it is found 0] that the critical current for switching is « 10 7 A/cm 2 , which is in very good 
agreement with experiments 0. 

The hysteresis loops for the case of the Co switching magnet consisting of a single atomic plane are 
reproduced in Fig. ^] The values of h p , 7, H u0 , and 9 are the same as in the previous example. However 
the ratio r « — 1 is now negative and the hysteresis loops in Fig. 1121 illustrate the interesting behaviour 
discussed in Section llTll when the system subjected to a bias higher than a critical bias moves to the "both 
unstable" region shown in Fig. [fjj As in Fig. [7|the points on the hysteresis loop in Fig. 1121 corresponding 
to the critical bias are labelled by asterisks. Fig. I12f b) and Fig. Ufa) are in close correspondence because 
Fig. Efa) is for r c < r < and in the present case r = —1, r c = —2/{^h v ) = —2. Also, from Fig. EJa), 
<?H < so that vg\\ in Fig. 13a) has the same sign as the voltage V in Fig. I12f b). 



VI. SUMMARY 



Spin-transfer torque is responsible for current-driven switching of magnetisation in magnetic layered 
structures. The simplest theoretical scheme for calculating spin-transfer torque is a generalised Landauer 
method and this is used in Section [H] to obtain analytical results for a simple model. The general 
phenomenological form of spin-transfer torque is deduced in Section IIIII and this is introduced into the 
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Landau-Lifshitz- Gilbert equation, together with torques due to anisotropy fields. This describes the 
motion of the magnetisation of the switching magnet and the stability of the steady states (constant 
current and stationary magnetisation direction) is studied under different experimental conditions, with 
and without external field. This leads to hysteretic and reversible behaviour in resistance versus bias 
(or current) plots in agreement with a wide range of experimental observations. In Section IIVI the 
general principles of a self-consistent treatment of spin-transfer torque are discussed and the Keldysh 
formalism for quantitative calculations is introduced. This approach to the non-equilibrium problem of 
electron transport uses Green functions which are very convenient to calculate for a realistic multiorbital 
tight-binding model of the layered-structure. In Section El quantitative calculations for Co/Cu/Co(lll) 
systems are presented which yield switching currents of the observed magnitude. 

This study of current-driven switching of magnetisation was carried out in collaboration with J. Mathon 
and A. Umerski and financial support was provided by the UK Engineering and Physical Research Council 
(EPSRC). 
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